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The angular momentum or vorticity of Cooper pairs is shown to affect strongly the charge transfer 
through a small superconducting (S) island of a single electron transistor. This interplay of charge 
and rotational degrees of freedom in a mesoscopic superconductor occurs through the effect of 
vorticity on the quantum mechanical spectrum of electron-hole excitations. The subgap quasiparticle 
levels in vortices can host an additional electron, thus, suppressing the so-called parity effect in the 
S island. We propose to measure this interaction between the quantized vorticity and electric charge 
via the charge pumping effect caused by alternating vortex entry and exit controlled by a periodic 
magnetic field. 

PACS numbers: 85.35.Gv, 73.23.Hk, 74.25.Ha 


There are two fundamental quantization phenomena 
which are manifested in different aspects of physics of 
superconducting (S) metals: (i) quantization of charge of 
superconducting carriers or Cooper pairs in units of dou¬ 
ble electron charge 2e and (ii) quantization of trapped 
magnetic flux in units of the flux quantum $o = he/ 2e 
which originates from the quantization of the angular 
momentum of Cooper pairs or vorticity. Basic supercon¬ 
ducting theory shows that these two quantization rules 
cannot be considered independently^ Certainly it would 
be exciting to clarify if this interplay of charge and vor¬ 
ticity quantization can reveal itself in possible magneto¬ 
electric phenomena inherent, e.g., to a superconducting 
state containing topological defects with quantized vor¬ 
ticity, namely, Abrikosov vortices. Unfortunately for the 
bulk samples it is extremely difficult to observe these ef¬ 
fects experimentally. The reason is that for typical metals 
the large value of the Fermi energy Ep almost completely 
suppresses all the electrostatic charge phenomena caused 
by vortices: the vortex core charge is small due to the 
very small ratio A /Ep ~ 10” 5 — 10 _2 j^^ Here A is the 
superconducting order parameter. 

In the present Letter we dare to make a suggestion 
overcoming the above difficulties based on the powerful 
methods provided by superconducting single electron de¬ 
vices (see, e.g., Refs. 0 and 0 and references therein). 
The key point of this idea is that by creating or remov¬ 
ing a single vortex in a small superconducting island, 
either the odd or even electron number will be favored. 
Thus the vorticity and the charge of the island will be 
coupled. To address electrons one by one we propose 
to use a single electron transistor (SET), i.e., a small 
metallic Coulomb-blockaded island with total electric ca- 



Figure 1. (color online) Setup of the NISIN SET with a bias 
voltage ±U/2 applied to the normal metal electrodes tunnel 
coupled to the central S disc with capacitances Cl and Cr. 
Magnetic field is applied perpendicular to the disc plane. 


pacitance C = Cp + Cr + C g (see Fig. |T]) coupled to 
the leads by tunnel contacts. Large Coulomb energy 
Ec = e 2 /(2C) of the island compared to the tempera¬ 
ture T prevents an additional electron to tunnel in and 
allows one to manipulate the charge state n of the island 
in a controllable way by varying the electrical potential 
of the gate electrode V g = en g /C g (see Fig. [I]), where C g 
is the gate electrode capacitance. For an S island this 
physical picture becomes more complicated due to the 
electron number parity effect^— This effect consists in 
the 2e - periodic dependence of the observables on an 
applied gate voltage and provides, thus, a direct confir¬ 
mation of the Cooper pair charge quantization. At finite 
temperatures the parity phenomena are controlled by the 
free energy difference SF = A — T\nN e ff between the 
states with odd and even number of excess electrons in 
the granule. Here N e ff equals the effective number of 
available states for an additional particle. It is clear that 
the parity effect is observable only for a positive value 
of this free energy barrier, i.e., at low enough tempera- 
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tures: T < T* = A/In N e ff. Applying an external mag¬ 
netic field one can suppress partially the S gap and, thus, 
suppress the parity phenomenon. One can observe this 
suppression either by measuring the change of periodicity 
of the SET characteristics vs the gate voltage in an ap¬ 
plied magnetic field IF,— or by using a varying magnetic 
field for controlled electron transfer through the SET at 
a fixed gate potential. At low temperatures the guaran¬ 
teed suppression of the parity effect can be achieved by 
introducing a vortex line in the granule which provides a 
natural trap for an entering electron. 



Applying an oscillating magnetic field, i.e., changing 
periodically the island vorticity, we can induce the even- 
odd transitions in the number n of electrons trapped on 
the island. Choosing the gate voltage as shown in Fig. [2] 
by a black dot, one can switch between the states n = 0 
and n = 1. Without the vortex, the even electron num¬ 
ber, n = 0, shown by the large white diamond is fa¬ 
vored. With the vortex, the odd electron number be¬ 
comes preferable as shown by the dashed red diamonds. 
Applying a constant bias voltage to the SET one can 
convert this modulation of the charge state into unidirec¬ 
tional charge pumping. The above picture of the vortex 
controlled parity effect can change at low temperatures 
less than the minigap wo < A in the spectrum of quasi¬ 
particles trapped in the vortex core. This new energy 
scale wo arises from the quantization of the spectrum of 
single particle excitations confined within the core by the 
Cooper pair potential?^ Though this minigap in most su¬ 
perconductors is small compared to the bulk gap it can 
still restore the parity effect. The free energy value paid 
for the even-odd transition in the electron number in the 
vortex state can be estimated as SF V = u >o — T hi N v , 
where N v ~ kpL , L is the length of the vortex line and 
fcp is the Fermi momentum. The condition SF V = 0 
gives us the temperature T* = u>o /In N v separating the 
regimes with e and 2e charge periods in the vortex state 
and, thus, at temperatures below T* the parity effect can 
be restored. It is quite useful to emphasize here a simple 
analogy with the parity effect in the Josephson junction^ 
where ujq should be replaced by the minigap that depends 
on the phase difference between the SC leads. 

We now proceed with the study of an exemplary SET 
setup (see Fig. [T|) which allows us to illustrate the above 
charge-vortex interplay. Hereafter we focus on the single 
electron transport between the normal metal leads and 
do not consider possible magnetic pumping based on the 
use of Cooper pair sluices in Josephson systems with S 
electrodes^ The size of the Coulomb blockaded S is¬ 
land is assumed to be of the order of several coherence 
lengths £ so that applying an external magnetic field we 
can introduce at least one vortex in this island. The elec¬ 
tronic transport through this device can be described by 
a standard rate equation accounting for parity effects? 14 
For the sake of simplicity we restrict to a two-level ap¬ 
proximation assuming low temperature regime T -C Ec 
and taking the gate voltage interval 0 < n g < 1. The 


Figure 2. (color online) (a) Stability diagram of the S island 
in the plane n g — V. Red dashed lines correspond to the 
normal state stability diagram. The zero current plateaus 
in the regime of the parity effect are shown by the white 
diamonds, (b) Ground-state energy of the S island vs gate 
voltage with (blue solid lines) and without (red dashed lines) 
parity effect. 

equation for the n = 1 charge state probability pi reads: 

= Fo-nPo -ri_>oPi , Po = l-Pi> (!) 
at 

where To-j-i and Ti^o are the rates for the electron tun¬ 
neling into and out of the island, respectively. These 
rates are, of course, determined by the sum of contri¬ 
butions coming from the transport through the contacts 
with left and right electrodes: r o-*-i = T e L [U L \ + T e R [U R ] 
andFi^o = T ° L [-U L \ + F° R [-U R \, where 

oo 

= / Vj(£)fN(s-U)[l-f§(e)]de (2) 

—OO 

is an increasing function of U. Here Rpj is the resis¬ 
tance of the jth tunneling junction, vj(e) is the local 
density of states (LDOS) of the island near the jth junc¬ 
tion normalized to its normal state value izjv( 0), index 
j = L,R stands for the left and the right junctions, 
/jv(e) = (e E / T + 1) _1 is the Fermi distribution function 
in the normal leads, and f§"' 1 (e) is the distribution func¬ 
tion in the S island describing the states with an even 
(odd) total number of electrons. The Coulomb blockade 
effect and the bias voltage Vl,r = ±V/2 determine the 
energy cost Ul,r = Ec(2n g — 1) — eV R ,R for tunneling. 

The increasing magnetic field and vortex entry affect 
both the LDOS vl,r{s) and distribution function fg°(e) 
in the above expressions. To find the distribution func¬ 
tion /g (o) (e) we assume that the zero (single) charge state 
corresponds to an even (odd) total number of electrons 
and use the so-called parity projection technique^ ' 15 ’ 16 

Mo), s _ Me) Texp(-2A 9P )/ a (£) 

/s U l±exp(-2 N qp ) ’ U 

where /F,s(e) = (e e ! 1 ± 1) _1 are the Fermi and Bose 
distribution functions. The number of quasiparticles can 
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be expressed as 


N qp = 2i//v(0) J dV J v(e, r)f F {e)de . (4) 

In the limit /f,b "C 1 one can neglect the difference be¬ 
tween these distribution functions and reduce © to the 
formii fs ’°( £ ) = 4l e,ofp { £ ) with the factor A e = A~ l = 
tanh(A r gp ). In the low temperature limit T <C £ m ; n with 
£ m in being the minigap in the quasiparticle spectrum of 
the island, we obtain N qp ss iV e //e _Emi,, / T , where N e ff 
is a slow function of temperature T (see Appendix [A] for 
details). 

Within the region of the essential parity effect (when 
\Ak —1| ~ 1) we can rewrite the tunneling rate as follows: 


Fj[U > 0] = Ij(U)/e 


rj P < o] = 


,U/T _ 1 


Ai. 


(5a) 

(5b) 


where the “seed” IV-characteristic of the tunnel junction 
in the absence of the Coulomb effects is 


OO 

W)= J ^[fN(e-U)-f N (e)]de . (6) 

— OO 


Note that Ij(—U) = —Ij(U). Further calculations should 
assume a certain model describing the dependence of the 
IV curves Ij(U) and the number of quasiparticles N qp 
on the applied magnetic field. For the sake of simplicity 
we consider the S island to be symmetric (see Fig. [TJ) 
assuming LDOS and tunnel resistances at both junctions 
to be equal, i.e., Vj = v(e) and Rtj = Rt- In this 
case key parameters governing the behavior of the IV 
curve, i.e., the minigaps e-j in the quasiparticle spectrum 
at the jth junction are also equal £j = £ e dge ■ The most 
important part of Ij(U ) controlling the charge transfer 
corresponds to small voltages (U < £ e dge) when the IV 
curve reveals the temperature activated behavior, 

OO 

Ij(U) ~ J- e -^ d9e ~ u )/T J v[s edge + Tx)e~ x dx . (7) 

o 

In the large voltage limit (U £ e d ge ) we assume a linear 
dependence Ij(U) = U/eRr ■ Note that we neglect here 
a low voltage contribution to the current arising from the 
exponential tail of the residual density of states localized 
inside the vortex core. 

Thus, the basic characteristics of our rate equation are 
determined by the magnetic field dependence of two en¬ 
ergy scales: (i) the spectral gap £ e d ge at the junctions 
and (ii) the minimal spectral gap £ m i n over the island. 
Considering an exemplary geometry shown in Fig. [l] one 
can see that the energy scale £ m i n is determined by the 
maximum of the local superfluid velocity vs reached ei¬ 
ther at the edge of the S disc or in the vortex core. The 



t/x 


Figure 3. (color online) The time-dependence of the energy 
scales £ e dge and e m in (a), the tunneling rates normalized to To 
(b), the probability distributions (c), and the instantaneous 
currents through jth junction (d). 


gap £ e d g e at the junctions is determined by the geometry 
of the S leads attached to the disc. Adding these S leads 
one can control the magnetic field effect on the tunnel¬ 
ing DOS and parity phenomenon independently. Taking, 
e.g., the diffusion limit with the coherence len gth well 
exceeding the mean free path ( we find (see Refs.ll£a-t20h: 

£edge = A(H)(l- 1 2 I l 3 ) 3 / 2 , ( 8 ) 

where A (H) = A(0)e -7r7ff / 4 , 7 h = ^( v s)/[ 2D A (if)], D 
is the diffusion coefficient, (u|) = (ttDHiv/^o) 2 /3, and 
w t; is the width of the S lead. Estimating now the 
energy scale £ m i n we can use the same expression © sub¬ 
stituting 7 jj = Hvg/[2DA(H)\ with the maximum local 
superfluid velocity vs- Assuming the screening effects to 
be small, i.e., when the disc radius is smaller than the 
effective London penetration depth, we get: max us = 
ttDHR/$ 0 . Considering clean limit £ <C i we should put 
£ m in = A (H)-hkFVs, where A (if) ~ A(0)[1— aH 2 / H 2 }. 
Here, H c ~ $ 0 /R£, is the field of the first vortex entry 
and a is a numerical factor of order unity. Just before the 
entry of the first vortex £ m i n tends to zero in the clean 
limit and remains at finite value in the dirty limit . 21-23 
After the vortex enters £ m in equals the minigap u >0 in the 
clean limit and turns to zero in the dirty regime. 

In order to model vortex induced pumping we assume 
the following time dependence of the spectral gaps with 
the period t = t 0 +t v (see Fig. © dictated by the piece- 
wise constant magnetic field applied: £ e dge = £min = 
A(0) for 0<t <t 0 and £edge(Rm) — A„, £min(-H m ) — £() 
in the interval to < t < t. The characteristic times of the 
vortex entry/exit are assumed to be negligible comparing 
to t v and to- Changing the value £0 from zero to wo we 
have the crossover from the dirty to the clean limit (at 
H m > H c ). The average current flowing through the j-th 
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Figure 4. (color online) (a) The current (J) averaged over 
the period vs V for H m = H c (green solid line) and IV curves 
at zero H = 0 (red dashed) and maximal H = H m (blue 
dasli-dotted) magnetic fields; (b) The current (J) vs squared 
amplitude of ac magnetic field H m normalized to the field of 
the first vortex entry H c at V = 0..0.5 Ec from left to right. 


junction at time instant t can be written as follows 

Jj(t) = e (T^UMPoit) - T°[-UMPi(t)) ( 9 ) 

using the solution of Eq. © 

Pi{t)=Pi(0)e-ti Tsdt '+ [ rWe-/*' r ^"dt', (10) 

J o 

Po(t) = 1 — Pi(t), where Ts = Eo-n + Ti^o- hi the case 
of periodic magnetic field protocol £ e dge(t+T) = £ e dge{t), 
e m i n (t + t) = e m in (t) the probability distribution arrives 
at the periodic steady solution after transient processes, 
when we can impose the condition pi(0) = Pi(t). As¬ 
suming Ur = Ul — \eV\ > 0 , Ul < A t , and Rl = Rr we 
can simplify the expressions for the tunneling rates m 


charge transmitted through both junctions should be 
equal to the electron charge. Due to the large ratio 
r^/r^j = e\ eV \/ T » 1 of the tunneling rates the most of 
this charge transfer e(l — exp(—|eV|/T)) occurs through 
the left junction, while only the exponentially small part 
of it ~ eexp(— |eVj/T) is transmitted through the right 
junction. The vortex exit should be accompanied by the 
discharge of the island which occurs with equal rates 
through the both junctions (Ulbl) . As a result, half of 
electron charge exits the island through each junction. 
Summing up the total charge transmitted through the 
system per cycle we find Q = e(l/2 — e~^ eV ^ T ). The 
above symmetry of the discharging processes results in 
a rather strong shot noise in the system: the fluctuat¬ 
ing transmitted charge equals to e /2 and the resulting 
current noise is given by the expression yj (<5 J 2 ) = ef / 2. 
The last two terms in m appear if the time intervals of 
two stages t v and to (with and without vortex, respec¬ 
tively) become comparable or shorter than the charac¬ 
teristic charging times l/T^, and l/T^- Therefore the 
maximum operation frequency / = 1 /r is limited by a 
single quasiparticle tunnel rate T 0 /N e ff. Besides the ef¬ 
fect of the frequency the average current (J) also deviates 
from ef /2 at small bias voltages and/or small magnetic 
field amplitudes, due to the dependence of total rates T|, 
on these parameters (see Fig.©. The terms proportional 
to 70 and originate from the leakages and lead to cur¬ 
rents exceeding e//2 shown at Fig. ©for larger V and/or 
H m values. Nevertheless satisfying the conditions 

A'e///r 0 < to < 7o 1 ) (13a) 

e (A v -U L )/Tj To < tv ^ 7 -l ) ( 13b ) 


r e 

1 L,R 

r° 

1 L,R 


P oe -( £ ed 9 e(i)-t/L,H)/T 

r 0 e- £ed ^)/ T 

tanh(A^ e j/e _£min ^)/ T ) ’ 


(Ha) 

(lib) 


where we neglect the slow time dependence of the pa¬ 
rameter To = e t ' Sedge( ' t ' l ~ u ^ T Ij(U > 0,t)/e (see Eq. ©). 
Substituting these expressions we find the current aver¬ 
aged over the period r = 1 //, 


(J) = - [ 1 - 2 e |ey|/T + yoto 


2 1 


ry;£o c — Fy^tv 


( 12 ) 


shown in Fig. ©a) and (b) vs bias voltage V and the mag¬ 
netic field amplitude H m in the time interval to < t < t. 
Here the total tunneling rates and in the vor¬ 
tex and Meissner states are determined mostly by the 
maximal rates F^ a: 2 Fo/lV e //, ~ Yoe ~ l ' A ' , ~ UL ' , / T , 

while 70 = r 0 e-( A -^)/ T , lv = are the 

small leakage rates. The first two terms in the above 
expression can be obtained from the following reason¬ 
ing. When the vortex enters the S island the total 


one can obtain the plateau of the average current at e //2 
(see Appendix [C] for details). These conditions can be 
met provided we set eV T, min(T, eo — TlnfV e //) <C 
Ul A — T\nN e ff. The breakdown of the above condi¬ 
tion at the lower bound of Ul can signal the presence of 
the minigap in the vortex core in the clean limit £q = ujo ■ 
Note that for non-symmetric case, the resulting average 
current will deviate from ef/ 2. However, a finite aver¬ 
aged current ranging between 0 and ef can be obtained 
in this case as well. 

To sum up, we have studied the interplay between vor- 
ticity and electric charge which can manifest in condi¬ 
tions of Coulomb blockade through the vortex induced 
suppression of the parity effect in mesoscopic samples. 
Vortex entry and exit from the sample is shown to be 
accompanied by synchronized entry and exit of a single 
electron charge. Applying the bias voltage and oscillating 
magnetic held one can observe a vortex governed turn¬ 
stile phenomenon: the switching between the Meissner 
and vortex states periodically opens the device for single 
charge transfer. Thus, we have demonstrated that the 
SET devices provide a unique tool for manipulating the 
collective dynamics of charge and vorticity in mesoscopic 
superconducting samples. 
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Within the region of the essential parity effect (when 
\Ak — 1| ~ 1) we can obviously neglect the term propor¬ 
tional to T qv and obtain the Eqs. © from the main text: 


Tj[U > 0] = Ij(U)/e 1 + 

■pfcrrr ^ nl Ij(U)/e a 
T J \y < - 1 _ e -U/T Ak 


Al 


,u/T _ 1 


(A7a) 

(A7b) 


Appendix B: Derivation of Eq. (1121) 


Appendix A: Derivation of expression for N qp and of 
Eqs. © 


In the low temperature limit T <C £ m ; n with £ m ; n being 
the minigap in the quasiparticle spectrum of the island 
Eq. (??) from main text can be written as follows 

N qp » N eff e~^ T , (Al) 


where 


N e ff = 2i//v(0)T J dV j u(£ min + T-x 1 r)e~ x dx (A2) 
is a slow function of temperature T. Indeed, 

dNeff 


dT 


= 2i/jv(0) J dV J xv(e m i n + T ■ x,r)e X dx~ 
Neff < i _ (A3) 


T T 


To derive Eqs. © from the main text it is convenient 
to rewrite Eq. © for the tunneling rates in the form: 


r i[ t/ ] = ^T / ^)f N {e-U)[l-f N {e)\de 

— OO 

oo 

+ J fN(e-U)-f N (e + U)} , 

3 0 

(A4) 

where we explicitly separate the parity effect contribu¬ 
tion. 

These expressions read 


r*[u] = 


1 - e~ u / T 


+ (Ak — 1) 


' IAW l)/e 

e |t/|/r _ i 


F e - |c/|/T 

qp g 


(A5) 


where Ij(U ) is the “seed” IV-characteristic of the tunnel 
junction in the absence of the Coulomb effects given by 
the Eq. (??) of the main text and 


F -_L 

qp ~ e 2 R 


o-ej/T 


Tj 


OO 

J v j ( e i 


+ Tx)e X dx . (A6) 


Using the magnetic field protocol considered in the 
main text £ e dge = £min = A(0) for 0 < t < to and 
£edge(N m ) — A y , £min(N m ) — £-0 hi the interval tfl <£ t <£ 
r, and the periodicity condition 


Pi(0) =Pi(t) 


(Bl) 


one can rewrite the solution (in in the main text as 
follows 


e (l — e 

Pl(t ^ to) Pl,eq (jPl,eq Pl,eq) ^ ^_pO £ Q _pp ^ 5 

" (B2) 

g — to) / 2 — g — 

Pi(t > t 0 ) = pl eq + (p° heq - p v heq ) - 1 _ e _ro tp _ r , tv —~ 

"(B3) 


Here the superscript 0 (v) corresponds to the time in¬ 
terval 0 < t < to (to < t < t, t = t 0 + tv) and 
Pi, eq = rg^i/rl is the adiabatic solution in the cor¬ 
responding time interval. 

Substituting this solution to Eqs. ©[13 from the main 
text one can obtain 


{J)t — Jeqto + Jeqty + e 




(Pl,eq - Pl,eq) 


(l - (l 


~r r.t v 


1 — e 




(B4) 


where + if^ and J k eq = e^T^o - 

is the current flowing through the island 
when the probabilities follow the adiabatic solution p k q . 

Using the expressions for all tunneling rates, assum¬ 
ing the conditions of Eqs. (fl3l) from the main text to be 
valid and keeping only the first order corrections in small 
parameters e ~\ eV ^ T , "foto, 7i Jv, e~ r x to , and e~ r x tv we 
obtain 


P 


o 

eq 



lv_ 7fc ~ p Tl 

pu > J eq ~ e r, ■ 


A Iv \eV\/T £R 

r s ~ 2r s ’ r° 


1 7o 

2 + 2r?,' 


(B5) 

(B6) 


Substituting these expressions into (IB4I) we get Eq. m 
from the main text. 


o 
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Appendix C: Range of parameters for current 
plateau observation 


In this section we consider ranges of the parameters 
where the plateau of the current Jj (t) averaged over the 
period r is close to ef / 2, i.e., (J) = e/(l ± e)/2 with a 
certain small £C 1. Using the result m from the main 
text one can roughly rewrite this condition as follows 


e-' ey ' /T ,7o to,7^,e- r ^,e- r ^ 


< 


(Cl) 


Further we focus on IV characteristic for (J)(V) as¬ 
suming that all other parameters (T, A, A„, eq, E c h = 
Eci^rig — 1 ), N e ff, to, and t v ) are chosen to be optimal 
for rather small bias voltages. 

Considering all the corrections to the current plateau 
one can separate them into two groups: (i) voltage- 
independent corrections 

2iy o/ iv.„ > indA), J h r ;;; / A ; , - T .(; T) <«, ^ 

and the necessary condition at V = 0 

r 0 e- (A - B ‘* ) / r t 0 < e, (C3) 


(ii) voltage-dependent corrections 


e~\ eV \/ T < e , Toe^ A - E ^ eV ^ T to < e, 

r 0 e~ (Av ~ Ech ~ leVl)/T t v > ln(l/e), (C4) 
Foe-^+^-I^D ! T t v < e, 


which can be rewritten as the conditions on the bias volt¬ 
age 


\eV\ > Tln(l/e) , (C5a) 

\eV\ > A V -E ch -T ln[r 0 t„/ ln(l/e)] , (C5b) 

\eV\ < A - E ch — T ln[r 0 f 0 /e] , (C5c) 

\eV\ <A V + E ch -Tln[T 0 t v /ln(l/e)] . (C5d) 


The last term in (IC4I) originates from the 
condition 7 v t v < e for the case Ur < 
-Tln[l/tanh(A e//e - E °/ T )] < 0. 

One can see that the increase in the minigap £q model¬ 
ing crossover between vortex minigaps in dirty and clean 
limits breaks first the last voltage-independent condition 
in Eq. (IC2I) . As a result, the plateau of the averaged cur¬ 
rent will be shifted to 7 v t v as a whole without change in 
the range of the bias voltage (IC5|) . 
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